Symmetry Analysis of Second Harmonic Generation at Surfaces of Antiferromagnets 
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Using group theory we classify the nonlinear magneto-optical response at low-index surfaces of 
fee antiferromagnets, such as NiO. Structures consisting of one atomic layer are discussed in detail. 
We find that optical second harmonic generation is sensitive to surface antiferromagnetism in many 
cases. We discuss the influence of a second type of magnetic atoms, and also of a possible oxygen 
sublattice distortion on the output signal. Finally, our symmetry analysis yields the possibility of 
antiferromagnetic surface domain imaging even in the presence of magnetic unit-cell doubling. 
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I. INTRODUCTION 



Optical Second Harmonic Generation (SHG) has been 
proven to be a very useful technique for the investigation 
of ferromagnetism at surfaces. The obvious question is 
if this technique can also yield some new information 
in the case of more general spin configurations, such as 
antiferromagnetic (AF) ordering. An experimental an- 
swer to this question has been provided by Fiebig et al. 
101 , who obtained a pronounced optical contrast from AF 
180° domains of rhombohedral bulk Cr203. The authors 
attributed this contrast to the interference of magnetic 
and electric dipole contributions, the latter being present 
only below the Neel temperature. Since it is known that, 
in cubic materials, within the electric dipole approxima- 
tion, optical SHG originates only from surfaces, inter- 
faces, or thin films, an important question is if SHG is 
also sensitive to antiferromagnetism at surfaces of cubic 
antiferromagnets. In this paper, we will show that the 
surface of a cubic material can lower the symmetry of an 
AF fee crystal (two-sublattice antiferromagnct) in a way 
similar to the trigonal distortion in a four sublattice an- 
tiferromagnct Cr203. Besides, even the imaging of AF 
domains is possible also for many cubic materials that 
exhibit unit-cell doubling. 

The first theoretical explanation of linear magneto- 
optic effects in ferromagnets has been given by Argyres 
1^ in the 50s. He used linear response theory for current- 
current correlation functions. His microscopic explana- 
tion was already based on the combination of spin-orbit 
and exchange coupling. Experimental techniques for the 
detection of AF domain walls using linear optics in some 
special geometries were elaborated a few years later 
The interior of the domains has been visualized in piezo- 
electric AF crystals using a linear magneto-optical effect 
Q. However, linear optical experiments suffer from mix- 
ing the desired signal with a contribution from other lin- 



ear effects, such as birefringence or dichroism. A review 
of linear optical experimental methods for the investiga- 
tion of AF domains is given by Dillon |^ . 

The observation of domain structure in antiferromag- 
nets is more complicated than in ferromagnetic materi- 
als since the reduction of the spatial symmetry is, unlike 
for ferromagnets, not linked to an imbalance in the oc- 
cupation of majority and minority spin states. On the 
basis of group theoretical considerations. Brown et al. 
^ proposed the use of linear optical effects, namely gy- 
rotropic birefringence, for the observation of AF domains 
related to each other by the space-inversion operation. A 
theoretical review of effects found by a group-theoretical 
approach is presented by Eremenko and Kharchenko [Q. 
They performed a comprehensive study of linear optical 
effects for various AF materials. Another effect proposed 
recently by Dzyaloshinskii et. al. ^ gives a possibility to 
detect antiferromagnetism taking advantage from optical 
path differences from antiferromagnetically coupled but 
intrinsically ferromagnetic planes. 

Nonlinear optics exhibits an additional degree of free- 
dom, since its elementary process involves three photons 
instead of two in linear optics. For that reason, some 
authors, e.g. Frohlich suggested the application of 
nonlinear optics even for k-selective spectroscopy, since 
multi-photon phenomena allow for the "scanning" of a 
small part of the Brillouin zone, at least for semicon- 
ductors. Recently, non-linear optics has attracted more 
and more attention for the investigation of magnetism 
due to its enhanced sensitivity to twodimensional ferro- 
magnetism | [Tc| ]. The magnetic effects are usually much 
stronger than in linear optics (rotations up to 90°, pro- 
nounced spin polarized quantum well state oscillations 
magnetic contrasts close to 100%) |l|,|l§|. An ex- 
ample of ferromagnetic effects measurable only by SHG 
deals with the existence of surface magnetism in very thin 
films of Fe/Cu(001) and is given in Ref. ||lq]. Nonlinear 
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optical effects were invoked to explain the behavior of 
lasers in magnetic fields |16|, t o investigate high temper- 
ature superconductors |l7yi^ ], and to study structures 
composed from alternately ferro- and antiferromagneti- 
cally ordered thin films One of the first theoret- 

ical investigations of the possibility to apply nonlinear 
optics to antiferromagnetism was performed by Kielich 
and Zawodny | pO[ . However, the first experiments con- 
cerning the detection of the AF domains in materials 
such as Cr203 were carried out only recently |^,^. Al- 
ready in the 70s, it has been proposed that experi- 
mental studies of dc magnetic and electric field-induced 
SHG could become an effective method of determining 
the crystal structure of solids, the symmetry of which 
cannot be investigated by other methods. Extending this 
idea towards surface crystallography provides us with a 
new technique for determining the spin configuration in 
a given surface structure. In turn, it permits to use a 
known magnetic configuration for the determination of 
the surface structure. All the mentioned effects are more 
difhcult or even impossible to obtain in linear optics, and 
moreover other linear methods like neutron scattering 
have difficulties to probe AF spin configurations. 

The nonlinear magneto-optical susceptibility tensor 
x'e^^ (the source for SHG within the electric dipole ap- 
proximation) has predominantly been investigated from 
the symmetry point of view. A classification following 
this approach, with tensors of a rank up to six, has 
been performed by Lyubchanskii et al. ||2^- ^ , |l3| , ^ . In 
Ref. the authors include the magnetization-gradient 
terms and apply the group-theoretical classification to 
higher-rank susceptibility tensors. This approach then 
allows them to study the thickness and the character 
(Bloch vs. Neel type) of domain walls. An attempt by 
Muthukumar et. al. to calculate the Xg;'^'' tensor el- 
ements for the antiferromagnetic Cr203 both from group 
theory as well as from the microscopic point of view is 
rather unique. They implemented a (Cr06)2 cluster, thus 
taking into account only half of the spins present in the 
elementary magnetic cell. In this approximation they 
explained the SHG from Cr203 as observed by Fiebig 
et al. ||f[ and they were able to give a quantitative esti- 
mate for that. Tanabe et al. | |29[ |, however, pointed out 
that the occurrence of purely real or imaginary values 
of the tensor elements plays a decisive role for the exis- 
tence of SHG from this substance. They found that for 
a (Cr06)2 cluster SHG can take place only in the case 
where the tensor elements are imaginary, and thus should 
vanish in Muthukumar 's approximation. They proposed 
to take into account the full unit cell with four inequiv- 
alent Cr ions including their "twisting" interaction with 
the environment. However Tanabe et al. neglected the 
dissipation in the process of SHG js^], which is a rather 
crude approximation. In general, taking into account the 
dissipation makes the xi?"'' tensor elements complex and 



invalidates their separation in purely real and imaginary 
ones | ]3l| . 

Lifting the inversion symmetry of a crystal is the source 
for SHG. Lyubchanskii et al. ]2^ , p6t suggested crystal 
lattice deformations and displacements as possible rea- 
sons for SHG from YIG films. In the case of Cr203 and 
YBa2Cu306-i-5, described by Lyubchanskii et al. p5| , p6[ , 
AF ordering lowers the symmetry of an otherwise cen- 
trosymmetric crystal. In this paper, however, we rely on 
the idea that, rather than lowering the crystal symmetry 
in the bulk, SHG may also result from the breaking of 
inversion symmetry at the surface of a bulk inversion- 
symmetric system. 

Magnetically active oxide layers are of importance for 
the construction of TMR (tunneling magnetoresistance) 
devices, where a trilayer structure is commonly used. 
The central layer of TMR devices consists of an oxide 
sandwiched between a soft and a hard magnetic layer 
(these two layers are often composed from the same ma- 
terial but of different thicknesses). For these technologi- 
cal applications it is necessary to develop a technique to 
study buried oxide interfaces. Such a technique can be 
SHG. One of the most promising materials for the men- 
tioned devices is NiO. However, to the best of our knowl- 
edge, the understanding of its detailed spin structure is 
scarce - even the spin orientation on the ferromagneti- 
cally ordered (111) surfaces is not known. The technique 
presented here can shed some light on that issue. 

Our paper is organized as follows: in Sec. II we present 
our methods for obtaining sets of nonvanishing xfi'^^ ten- 
sor elements. In Sec. HI we present the results of our 
analysis, first for the nondistorted surface of a simple fee 
structure (subsection III. A), then for the the distorted 
one (III.B). Subsequently, we discuss the influence of a 
second kind of magnetic atoms (III.C) and of oxygen sub- 
lattice distortion (III.D). The issue of domain imaging 
is addressed in subsection III.E. Possible experimental 
geometries allowing for the detection of the mentioned 
structures and effects are discussed in Sec. IV. The con- 
clusions are presented in Sec. V. 

II. THEORY 

Based on group theory, Dahn et al. proposed a new 
nonlinear magneto-optic Kerr effect (NOLIMOKE) at 
the surface of cubic antiferromagnets. They also gave an 
example of an antiferromagnetic structure (NiO) and an 
optical configuration, where this new effect could be ob- 
served. Here, we perform a complete group-theory based 
analysis of coUinear AF fee low-index crystal surfaces. 
Surfaces of other crystal structures are as well described 
by our theory provided they are similar to fee crystal sur- 
faces, i.e. squares or hexagons. The results can be used to 
detect the magnetic order of a specific surface under in- 
vestigation and allow for the determination of the surface 
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spin configuration in some important cases. However, in 
order to calculate the SHG yield quantitatively, it is nec- 
essary to go beyond the present study and use electronic 
calculations of the nonlinear susceptibility. Group the- 
ory can give a unified picture of different experimental 
observations and predict new effects while the mi- 
croscopic origins of the observed phenomena may remain 
unclear. In order to be clear with respect to the essential 
notion of time reversal we would like to emphasize the 
point of view taken in this paper in the beginning. Here, 
we do not divide xi;'^'* into even and odd parts in the 
magnetic order parameter. Instead, the behavior of x[?'^'* 
with respect to the magnetic order parameter (which for 
ferromagnetic materials corresponds to the dependence 
of Xd,^"^ '-'^ magnetization) is fully taken into account by 
the considerations of the magnetic point group. At no 
stage of our consideration we invoke the notion of time 
reversal, consequently we do not apply the characteriza- 
tion of the susceptibility as c-tensor (changing its 
sign in the time-reversal operation) or i-tensor (invariant 
under the time-reversal operation) |]3l|] . 

Before we start our group theoretical classification of 
the nonlinear optical susceptibilities of AF surfaces we 
would like to emphasize the following four important 
points: 

(i) We are not interested in effects resulting from the op- 
tical path difference from adjacent crystal planes which 
are ferromagnetically ordered but only antifcrromagnct- 
ically coupled to each other. We do not consider this as 
an intrinsic AF effect. 

(ii) Cubic crystals that we are interested in reveal a center 
of inversion in the para-, ferro-, and all antiferromagnetic 
phases. Thus, within the electric dipole approximation, 
the SHG signal from the bulk vanishes. 

(iii) While in principle linear optical methods can be sen- 
sitive to the presence of a spin structure, in practice they 
are not useful because, within the group theoretical ap- 
proach, they cannot distinguish the AF phase from either 
paramagnetic or ferromagnetic, nor can they distinguish 
different AF configurations from each other. They have 
to resort to methods like lineshape analysis, where no 
strong statements characteristic for symmetry analysis 
can be made. 

(iv) Although the tensor elements for all the magnetic 
point groups are known and tabulated in the literature 
(e.g. Q), the connection between the different spin con- 
figurations described by us and the mentioned symme- 
try groups has not been made, except for some easy 
cases [Q. Thus, for SHG from antiferromagnetic sur- 
faces there has been up to now no connection between 
the group theoretical classification and the real situations 
found in experiments. 

The following part of the text should explain the fun- 
damentals of applying NOLIMOKE observations to in- 
vestigate antiferromagnetism of surfaces. 



Now we turn to SHG, the source of which is the non- 



linear electrical polarization P^f^' given by: 



(1) 



Here, E'^'^^ is the electric field of the incident light, while 
Xgj'^'' denotes the nonlinear susceptibility within the elec- 
tric dipole approximation, and eo is the vacuum permit- 
tivity. The intensity of the outgoing SHG light is l35| : 
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where Iq is the intensity of the incident light, F{) (/()) 
describe Fresnel and geometrical factors for the incident 
(reflected) light, i) and O angles of incidence and reflec- 
tion, respectively (??=0), and {ip) is output (input) 
polarization angle. According to Neumann's principle, 
"any type of symmetry which is exhibited by the crystal 
is possessed by every physical property of the crystal" 
P^ . To examine these physical properties, we determine 
the magnetic point group of the crystal lattice, thus de- 
termine its symmetries. The same symmetries must leave 
the investigated property tensor (in our case the nonlin- 
ear electric susceptibility X^.^^) invariant. This fact is 
mathematically expressed by the following condition: 

(3) 

Here, Z„,„'(n = i,j,k,n' = i',j',k',) is a representation 
of an element of the magnetic point group describing the 
crystal. For symmetry operations including the time re- 
versal there should be an additional "±" sign in Eq.(||), 
but we do not use it here since we exclude the time rever- 
sal from our consideration. In particular, from Eq.(^) it 
follows immediately that polar tensors of odd rank (such 
as x^ei^^ ) vanish in inversion symmetric structures. This 
explains why SHG is possible only at surfaces and in- 
terfaces, where this symmetry is broken. For a given 
spin configuration we apply Eq. (^ for every symmetry 
operation exhibited by the system. Thus, each of these 
symmetries gives rise to a set of 27 equations with 27 



unknown elements of the tensor Xei 
reduced to 18 equations, since 
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This set can be 



(4) 



which expresses the equivalence of the incident photons 
of frequency to, see also the reduced notation in Eq. (2). 
The analytic solution of even this reduced set of equa- 
tions seems cumbersome, but the set can be split into 
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several decoupled subsets. For example, an obvious sub- 
set in every case is the equation Xzzz = Xzzz-, this tensor 
element occurs nowhere else. The rank of other subsets 
is, for our cases, never higher than six. In this manner, 
one may obtain a set of forbidden elements of the suscep- 
tibility tensor as well as relations between existing ones. 

III. RESULTS 

First, we will define the notions of "phase", "case", 
and "configuration", used henceforth to classify our re- 
sults. "Phase" describes the magnetic phase of the mate- 
rial, i.e. paramagnetic, ferromagnetic, or AF. Secondly, 
the word "configuration" is reserved for the description 
of the magnetic ordering of the surface. It describes 
various possibilities of the spin ordering, which are dif- 
ferent in the sense of topology. We describe up to 18 
AF configurations, denoted by little letters a) to r), as 
well as several ferromagnetic configurations, denoted as 
"ferrol", "ferro2", etc. The number of possible configu- 
rations varies depending on surface orientation. Thirdly, 
we describe different "cases", i.e. additional structural 
features superimposed on the symmetry analysis. "Case 
A" does not have such additional features. In "case B" 
we address distortions of the lattice. "Case C" deals with 
two kinds of magnetic atoms in an undistorted lattice. 
In "case D" we take into account a distorted sublattice 
of nonmagnetic atoms, keeping the magnetic sublattice 
undistorted. All the analysis concerns coUinear antifer- 
romagnets, with one easy axis. 

The tables show the SHG response types for each con- 
figuration. The various response types are encoded by 
a "key", which is then decoded in Tab. |. This table 
presents the symmetries, domain operations, and non- 
vanishing tensor elements for each response type. This 
is done in order to shorten the overall length of tables, 
because a given response type can appear in several dif- 
ferent cases. 

Several spin structures depicted in Fig. |] and Fig. ^ 
are distinct configurations only in case B, and they are 
addressed in the tables that concern only this case. For 
the rest of the cases they are domains of other, fully 
described configurations, thus they are left out in these 
cases. The philosophy of the paper is that, to save some 
space, we show the spin structure in one figure for each 
surface (Fig. |^, IJ, and ||) for all the four cases (A-D), 
and depict the effects taken into account in the cases 
B-D only for the paramagnetic phase (Fig ^, 0, and ||). 
Table | also contains the information on the parity of the 
nonvanishing tensor elements: the odd ones are printed 
in boldface. In some situations an even tensor element 
(shown in lightface) is equal to an odd element (shown 
in boldface) , this means that this pair of tensor elements 
is equal in the domain which is depicted on the corre- 
sponding figure, but they are of opposite sign in the other 



domain. This happens in the structures where two pairs 
of domains are possible (two distinct entries in Table |). 
The tensor elements that change their parity in the do- 
main operation which is the inverse of the displayed one 
are shown in italic font. For example, the entry j) of 
Table | shows a tensor element xxx^ which is even under 
the operation 4^, this means that this tensor element is 
odd under — 4^. This strange at the first sight behavior of 
tensor elements is caused by the fact that under these op- 
erations, tensor elements are not mapped on themselves. 
In our example, after applying 4^ the tensor element xxx 
becomes yyy, without changing its sign. If we now ap- 
ply — 4^, yyy (which is now even) becomes xxx^ again 
without changing the sign. 

The parity of the elements has been checked in the op- 
erations 2^, 4z, and in the operation connecting mirror- 
domains to each other (for the definition of the mirror- 
domain structure see subsection E). The domain opera- 
tion(s) on which the parity depends is (are), if applica- 
ble, also displayed in this Table. If two or more domain 
operations have the same effect, we display all of them 
together. To make the Table || shorter and more easily 
readable some domain operations (and the correspond- 
ing parity information for the tensor elements) are not 
displayed, namely those that can be created by a super- 
position of the displayed domain operations. We also do 
not address the parity of tensor elements in the 6z nor 
3z operations for (111) surfaces nor any other operation 
that "splits" tensor elements, although these operations 
also lead to a domain structure . As will be discussed 
later (subsection E) it is possible to define a parity of the 
tensor elements for the 3^ and 6^ operations, however the 
tensor elements then undergo more complicated changes. 
The situations where the parity of the tensor elements 
is too complicated to be displayed in the Table are in- 
dicated by a hyphen in the column "domain operation" . 
For some configurations, none of the operations leads to 
a domain structure - in those configurations we display 
the information "one domain" . The reader is referred to 
the Appendix for the particularities of the parity check. 

As far as the first layer is concerned, we address all the 
spin configurations of the low index surfaces of fee anti- 
ferromagnets, with magnetic order vector lying in plane 
or perpendicular to it and antiferromagnetic coupling 
between nearest neighbors. For the (001) surfaces we 
also discuss the configurations, where the antiferromag- 
netic coupling exists between the second-nearest neigh- 
bors (configurations a), b), c), f), and o), along with d), 
g), and h) for case B.). We do not consider the coupling 
to the third and further neighbors. This would not give 
rise to configurations of different symmetries in two di- 
mensions. It may at most replace spins by grains (blocks) 
of spins in the configurations described by us. 

Throughout this paper we take into account the spin 
structure only of the first (uppermost) atomic layer. This 
is sufficient to study all the symmetries of (001) and 
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(110) surfaces both in the paramagnetic and ferromag- 
netic phases. For the (111) surface it is necessary to 
recognize the atomic positions (but not the spins) in the 
second layer for the same purpose. For the sake of com- 
pleteness we also present a study of (111) surfaces with- 
out this extension. However, in the antiferromagnetic 
phase, the spin structure of the second and deeper layers 
plays a role in determining the symmetry of the surface. 
This is presented in this paper using the (001) surface 
as an example. For the (110) and (111) surfaces it will 
be published elsewhere These structures can serve 
as simple models for deriving predictions for more com- 
plicated cases, while the full consideration of the second 
layer would not bring any new interesting results. Tak- 
ing into account the spin structure of the second layer 
(deeper layers do not bring up anything new to the anal- 
ysis) results in creating several (up to two for the (001) 
surface and three for the (111) surface) configurations 
out of each one addressed here by us. The symmetry of 
these configurations may remain the same or be lowered 
(sometimes even below the symmetry of the ferromag- 
netic phase) with respect to the "two-dimensional" con- 
figurations they are generated from. Consequently the 
distinction of the configurations from each other may be 
limited, but the possibility to detect the magnetic phase 
is not severely affected. Also our remarks on domain 
imaging remain valid. However the number of domains 
is increased, thus the possibility to identify each of them 
might be hampered. 

Consequently, one can state that the symmetry of an 
AF surface depends on two atomic layers. They are also 
necessary (and sufficient) to define AF bulk domains. As 
will be presented in our results, SHG can probe both 
these layers on AF surfaces. 



A. Equivalent atoms 

The predicted new nonlinear magneto-optical effects re- 
sult from the fact that the magnetic point groups of 
antiferromagnetic configurations are different from those 
describing paramagnetic or ferromagnetic phases of the 
same surface. Since, depending on the magnetic phase, 
different tensor elements vanish, it is possible to detect 
antiferromagnetism optically by varying the polarization 
of the incoming light. 

The current subsection discusses nonvanishing ele- 
ments of the nonlinear susceptibility tensor for an fee 
crystal consisting of only one kind of magnetic atoms. 
The influence of nonmagnetic atoms in the material will 
be discussed later. The configurations considered here 
are "ferrol", "ferro2", "ferro4", a), b), c), e), f), i), k), 
m), o), p), and r) for the (001) surface (see Fig. |^), 
"ferrol", ferroS", "ferro5", a), c), f), i), and k) for the 
(111) surface (see Fig. ||), and all configurations depicted 



in Fig. ^ for the (110) surface. Other depicted spin struc- 
tures form domains of these configurations and are not 
referred to in this subsection nor in the tables concerning 
the current subsection psf . 

All possible configurations (confs.) of a fee (001) sur- 
face are shown in Fig. [T], which displays the conventional 
rather than magnetic unit cells. However, these are suf- 
ficient to fix the spin configuration of the whole surface 
imposing of the following "convention" : the fee surface is 
constructed from the depicted plaquette in the way that 
neighboring spins along the x and y directions point the 
same way (alternate) if they are parallel (antiparallel) 
on the plaquette in these two directions. The spins in 
rows and columns where only one spin is presented are 
continued in the same way as the corner spins. For in- 
stance in the configuration a) of the (001) surface, both 
the right-hand side and left-hand side neighbors of the 
"central" spin will point upwards, while the spin direc- 
tion will be alternated along the x axis. This convention 
will be maintained henceforth (for a (111) surface one 
has to alter or keep the spins along three axes, instead of 
two). The smallest set that gives a complete idea about 
the spin structure is presented in Fig. || this "mag- 
netic primitive cell" does not give a clear picture of the 
crystal symmetries, however. The whole crystal lattice 
can be reproduced by translations of this cell, without 
performing other operations such as reflections or rota- 
tions. 

The SHG response types for the (001) monolayer are 
given m Table ||, for the paramagnetic, ferromagnetic, 
and all AF phases. We can observe several sets of allowed 
tensor elements. The Conf. r) will produce the same 
signal as the paramagnetic phase. The Conf. "ferrol" 
reveals a completely different, distinguishable set of ten- 
sor elements. In addition, the conf. "ferro2" produces 
another set of tensor elements, different from any other 
configuration. It is equivalent to the conf. "ferrol" ro- 
tated by 45°. In the confs. a), b), e), and o) we find 
the same tensor elements as for the paramagnetic phase. 
However, due to the lower symmetry, their values are no 
longer related to each other. Confs. c) and f) bring new 
tensor elements, thus allowing for the distinction of these 
confs. from the previous ones. Confs. i), k), m), p) reveal 
the same tensor elements as c) and f) but some of these 
elements are related. Thus one may possibly distinguish 
these two sets of configurations. Conf. "ferro4" presents 
a completely different, distinguishable set of the nonva- 
nishing tensor elements. Consequently, in six configura- 
tions (i.e. c), f), i), k), m), and p)) some susceptibility 
tensor elements appear only in the AF phase, allowing 
for the detection of this phase by varying the incident 
light polarization, as will be outlined in Sec. IV. In ad- 
dition, all other antiferromagnetic configurations but r) 
reveal the breakdown of some of the relations between 
the different tensor elements, compared to the paramag- 
netic phase, and thus can be detected as well. Gener- 
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ally, all the phases can be distinguished from each other. 
There exists as well a possibility to distinguish different 
AF configurations provided the corresponding tensor el- 
ements can be singled out by the proper choice of the 
experimental geometry. 

For the sake of completeness, we now present a short 
study of the (001) surface where the spin structure of the 
two topmost atomic layers is taken into account. The 
paramagnetic phase and all the ferromagnetic configu- 
rations remain unchanged with respect to the results of 
the previous paragraph (for the (001) monolayer). How- 
ever, most of the AF configurations previously addressed 
break up into two different configurations (sometimes 
even with a different symmetry). These configurations 
are constructed from the ones of the previous paragraph 
by assuming that the structure of the second atomic layer 
is identical with that of the topmost one but shifted along 
the positive x axis (indicated by x after the name of 
the original configuration) or positive y axis (indicated 
by y after the name of the "parent" configuration) in a 
proper way to form a fee structure; if only one configu- 
ration can be produced in this way we use the name of 
the original one. This construction is depicted in Fig. 
^, along with the corresponding conventional unit cells 
for the two topmost layers of the AF fee (001) surface. 
The resulting SHG response types are presented in Ta- 
bic III, In general, seven types of response are possible. 



Firstly, the paramagnetic phase reveals a characteristic 
set of tensor elements. Thus it can be unambiguously 
distinguished from any other magnetic phase. Secondly, 
confs. "ferrol", ax), ox), bx), by), ex), and ey) bring 
some additional tensor elements into play. The symmetry 
of confs. ax) and ox) is slightly different from the one of 
the rest of this group, since the mirror plane is rotated by 
90° around the z axis. A different set of tensor elements 
is brought up by confs. "ferro2" , i), m), and p). The 
difference between the response yielded by conf. i) and 
the other confs. in this group, due to a slightly different 
symmetry, can be compensated by rotating the sample 
by 90° around the z axis. Another, characteristic set of 
tensor elements is presented by conf. "ferro4" alone. The 
fifth type of SHG response is given by confs. ay), oy), 
and r). Tensor elements, that do not vanish in these con- 
figurations, are the same as for the paramagnetic phase 
but some relations between them are broken due to a 
lower symmetry in the AF phase. Confs. cx), fx), and 
fy) yield all tensor elements in an unrelated way. The 
last, characteristic type of response is presented by conf. 
k) alone. Consequently, the detection possibilities of an 
antiferromagnetic bilayer are slightly worse than those 
for a monolayer. Especially, a difficulty in distinguishing 
the ferromagnetic phase from the antiferromagnetic one 
may arise for some configurations where then the com- 
bination of SHG with linear magneto-optics is definitly 
required. There exists a possibility to distinguish AF 
configurations from each other, similarly to the previ- 



ous situation. In most configurations, the difference (in 
terms of the SHG response) between the bilayer struc- 
ture described here and the previously addressed (001) 
monolayer can be detected. 

We now turn to the (110) surface (Fig. which, in the 
paramagnetic phase, reveals a lower symmetry than the 
(001) surface. On the other hand, the number of sym- 
metry operations in the AF configurations is comparable 



to the (001) surface. In addition, as shown in Table IV, 
the resultin SHG response types are not very characteris- 
tic, so the detection possibilities for this surface are very 
limited. In particular, confs. a), b), c), g), h), i), j), k), 
and 1) give the same tensor elements as the paramagnetic 
phase. Confs. d), e), f), and "ferro3" bring new tensor 
elements. Other ferromagnetic configurations ("ferrol" 
and "ferro2") present different sets of new tensor ele- 
ments, making these configurations distinguishable from 
the others as well as from each other. Conf. "ferro4" 
yields a completely different set of tensor elements, how- 
ever this set is related to the one of conf. "ferrol" by 90° 
rotation. 

The study of the (111) surface (see Fig. ^ has to be 
separated in two subcases, according to whether we take 
into account only one atomic monolayer or more. In both 
subcases, we consider the same configurations. The SHG 
response types for the first subcase are listed in Table 0, 



and for the second subcase in Table VI. For the first 



subcase, confs. a), i), and k) reveal the same tensor el- 
ements as the paramagnetic phase, however due to the 
lower symmetry their values are not related to each other. 
Configurations c) and f ) present new tensor elements. As 
for the previous surfaces, the ferromagnetic phase reveals 
completely different sets of tensor elements, and the three 
ferromagnetic configurations can be distinguished from 
each other since they bring different tensor elements into 
play. Unlike for the (110) surface, the axes x and y are 
not topologically equivalent, and thus the fact that ten- 
sor elements of "ferrol" are related to those of "ferroS" 
by 90° rotation does not affect the possibility to distin- 
guish these two configurations. The ferromagnetic conf. 
"ferroS" brings up the same tensor elements as AF confs. 
c) and f), but the relations between the elements are 
different. The second subcase (more layers taken into 
account) gives different sets of allowed tensor elements 
(compared to the first subcase) for each but the "ferro3" 
configuration. Confs. a), i), k), and "ferroS" share the 
same set of allowed tensor elements and can be easily 
distinguished from the paramagnetic phase. Confs. c), 
f), and "ferrol" reveal all tensor elements, with their val- 
ues unrelated. Similarly, conf. "ferroS" presents another, 
distinguishable set of tensor elements. The possibility to 
distinguish the magnetic phases is rather limited. 

The symmetry analysis of nonvanishing tensor ele- 
ments for ferromagnetic surfaces in the case A have been 
performed by Pan et. al. ]To| |. Our analysis yields the 
same results, taking into account the corrections made 
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by Hiibner and Bennemann |4C 



B. Distortions of monoatomic lattice 

The rhombohedral distortion of the atomic lattice, de- 
scribed here and shown in Fig. makes the x and y 
axes of the (001) surface inequivalent, even in the para- 
magnetic phase. On the (111) surface, the y axis is not 
equivalent any longer to other axes connecting the near- 
est neighbors, namely 5'(xi/) and S'j^^.j,) (for the defini 



tion of the "S" and "H" axes see Fig. |5| the paramag- 
netic conf.). These inequivalences of axes are the rea- 
sons for the reduction of the number of symmetry op- 
erations in the paramagnetic phase. Because of this re- 
duction some spin structures that previously formed dif- 
ferent domains of a single configuration now cannot be 
transformed into each other and become "independent" 
configurations. This happens for almost every of the pre- 
viously addressed configurations of the (001) and (111) 
surfaces. Consequently, all the depicted spin structures 
are in fact configurations, and are addressed in this sub- 
section. 

The resulting SHG r espo nse types for the (001) sur- 
face are listed in Table VII. For this surface, only two 



of the ferromagnetic configurations, namely "ferrol" and 
"ferro2" can be easily distinguished from both the para- 
magnetic as well as the antiferromagnetic phases. These 
ferromagnetic configurations can be also distinguished 
from each other. On the contrary, all the AF configu- 
rations yield only two types of response, and in addition 
one of them is equivalent to the response of the paramag- 
netic phase. Consequently, it will not be possible to de- 
termine the surface spin structure, and the distinction of 
the AF phase from the paramagnetic one can be success- 
fully performed only in confs. a)-h) and o). Compared 
to the case A, there is an important symmetry breaking 
for most configurations. Thus, the distinction between 
the two cases (A and B) is possible (compare Tabs. || 
and[vi|). 

All the (110) surfaces of an fee crystal with a rhom- 
bohedral distortion are topographically equivalent to the 
(110) surface of the case A. The distortion only stretches 
the X or y axis, so the structure remains rectangular. 

The analysis of the (111) surface (depicted in Fig. |^) 
in the subcase of only one monolayer reveals sets of sym- 
metries very similar to the (110) surface, as it follows 
In fact, the (111) surface of a fee 



from the Table VIII 



crystal with a rhombohedral distortion can be treated 
as two rectangular lattices superimposed on each other. 
In turn, due to the distortion, it is not convenient any 
longer to describe the spin structures using "S" and "H" 
axes. The possibility to distinguish AF configurations is 
very poor, and two of the AF configurations (a) and k)) 
yield the same signal as the paramagnetic surface. In 



confs. b) - j), 1), and m) the AF phase can be distin- 
guished from the paramagnetic one, but they give the 
same signal as conf "ferro5" . Conf. "ferro2" can be eas- 
ily distinguished since it reveals a characteristic set of 
(all) tensor elements. Confs. "ferrol" and "ferroS" yield 
different sets of tensor elements, but they are related to 
each other by 90° rotation. Most of the configurations 
allow for the distinction of the cases A and B (compare 
Tabs. and |vn| . 

In the subcase of two monolayers of the (111) surface, 
the symmetry is dramatically reduced (see Tab. IX). 



Even in the paramagnetic phase the group of symme- 
tries consists of only one nontrivial operation, and this 
appears to occur also in the AF configurations a), i), k), 
and "ferro3". In all the other configurations all tensor 
elements are allowed due to the lack of any symmetry. 
Only confs. paramagnetic and "ferro5" allow for the un- 
ambiguous distinction of the cases A and B (compare 
Tabs. VI and IX). Consequently, this surface is not very 
useful to an analysis of the magnetic structure, with the 
exception of stating the distortion itself. 

As the conclusion of the case of the distorted sublat- 
tice of magnetic atoms, the surfaces give extremely lim- 
ited possibilities to investigate the magnetic properties. 
In our further study, we will limit ourselves to lattices of 
undistorted magnetic atoms. 



C. Structure with nonequivalent magnetic atoms 

We assume now that not all the magnetic atoms in the 
cell are equivalent. An example of such a structure is 
a material composed of two magnetic elements, but also 
a situation when the magnetic lattice sites are inequiva- 
lent due to different bonds to a nonmagnetic sublattice; 
distortions of the sublattice of nonmagnetic atoms that 
preserve the center of twodimensional inversion produce 
the same effect. Other distortions of the sublattice of 
nonmagnetic atoms will be discussed in subsection D. 
The magnetic moment at the distinguished positions can 
be changed or not - this does not affect the results ob- 
tained by symmetry analysis. The configurations consid- 
ered here are "ferrol", "ferro2", "ferro4", a), b), c), e), 
f), i), k), m), o), p), and r) for the (001) surface (see 
Fig. |), "ferrol", ferro3", "ferroS", a), c), f), i), and k) 
for the (111) surface (see Fig. ||), and all configurations 
depicted in Fig. ^ for the (110) surface. Other depicted 
spin structures form domains of these configurations and 
are not referred to in this subsection nor in the tables 
concerning the current subsection. 

The structure is depicted in Fig. ^ For the sake of 
brevity, we show the structure of the distinguished atoms 
only for the paramagnetic phase. All the configurations 
are the same as in case A, for all surface orientations. 
The already mentioned "convention" of alternating (or 
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not) spin directions along certain axes is applied regard- 
less of the atom type. This allows us to obtain the whole 
crystal surface from the small displayed fragment. 

Our analysis starts with the (001) surface of an fee 
crystal. The SHG response types for each configuration 
are listed in Table |^. In general, we can observe seven 
types of response. The first of them is represented by the 
paramagnetic phase alone. The second type of response, 
exhibited by the ferromagnetic "ferrol" and the AF a), 
b), e), o) confs., differs from any other type by some ten- 
sor elements. The confs. a) and o) reveal different tensor 
elements than the other configurations from the men- 
tioned group. However, the signal from confs. a) and o) 
is the same as for the confs. b), e), and "ferrol" if one 
exchanges the axes x and y. Thus, if the directions of the 
spins cannot be determined by another method, confs. a) 
and o) cannot be distinguished from b), e), and "ferrol". 
The next type consists of conf. f) and reveals all tensor 
elements, while no relations between them are enforced 
by the symmetry analysis. A completely different type 
of response is presented by conf. c) alone. Another type, 
where confs. i), m) and p) belong to brings the same ten- 
sor elements as conf. c), but there exist more relations 
between the elements due to a higher symmetry in these 
configurations. The next type is given by confs. "ferro2" 
and k) . As in conf. f ) all the tensor elements are present 
but this time there are some relations between them. In 
addition, confs. r) and "ferro4" yield a completely new 
set of tensor elements due to the preserved fourfold rota- 
tional symmetry. 

Thus, assuming one atom as distinguished may reduce 
the symmetry. New tensor elements appear in confs. a), 
b), e), f), k), o), and r) compared to case A (compare 
Tabs. H and In these configurations it is there- 

fore possible to distinguish the cases of equivalent and 
nonequivalent magnetic atoms, provided the tensor ele- 
ments that make the cases different can by singled out 
by the experimental geometry. There exists also a pos- 
sibility to distinguish different AF configurations in case 
C. The antiferromagnetic phase can be undoubtelly de- 
tected in the surface configurations c), f), i), m), and 
P). 

For the (110) surface, there are more possibilities to 
distinguish the configurations with nonequivalent mag- 
netic atoms than in the case A. However, the configu- 
rations still produce ambiguous signals (see Tab. |x| ). 
Confs. b), c), h), i), k), and 1) are equivalent to the para- 
magnetic phase. Conf. a) is equivalent to the ferromag- 
netic "ferrol" configuration, and conf. d) to "ferro2" . In 
addition, the confs. e), f), and g) are equivalent to the 
conf. "ferroS" and conf. j) gives the same signal as conf. 
"ferro4" . Even the presence of nonequivalent atomic sites 
in the lattice cannot be detected by SHG on this surface, 
since the symmetry of the (110) surface is usually not 
lowered further by the existence of equivalent magnetic 
sites (compare Tables |^ and X] . The only exception are 



the confs. a), d), g), and j) which give different tensor 
elements in the two cases. As in the case of equivalent 
atoms, the (110) surface is not very useful for the analy- 
sis. 

The study of the (111) surface must again be divided 
in the two subcases of one or more monolayers, respec- 
tively. Fig. 1^ depicts the situation in the paramagnetic 
phase. The SHG response types are listed in Tables XII 



and Kill for the first and the second subcase respectively. 



In the first subcase (one monolayer) the symmetry es- 
tablishes six different types of nonlinear response. The 
"paramagnetic" type (for the paramagnetic configuration 
only) is characteristic - all the other configurations have 
additional tensor elements. The next type of response 
(the ferromagnetic conf. "ferrol" and the antiferromag- 
netic conf. a)) brings some new tensor elements. Other 
tensor elements appear in the conf. k). Configurations 
"ferroS" and i) show another set of nonvanishing tensor 
elements. The confs. c) and f) reveal all tensor elements 
in an unrelated way. In addition, conf. "ferro5" presents 
a characteristic set of tensor elements. 

In the second subcase, only four different SHG re- 
sponses are possible. Firstly, the paramagnetic phase 
is characteristic - all the other configurations bring ad- 
ditional tensor elements into play. The next type of re- 
sponse is presented by confs. "ferroS" and i) - they yield 
some additional tensor elements. Confs. "ferrol", a), c), 
f), and k) reveal all tensor elements and no relations be- 
tween them appear from our symmetry analysis. Again, 
the conf. "ferro5" presents a unique set of nonvanishing 
tensor elements. 

Consequently, for the (111) surface, the symmetry 
breaking due to the presence of a second kind of mag- 
netic atoms has even more important consequences than 
for the (001) surface. In the situation of only one mono- 
layer, the distinction between the cases may be po ssibl e 
for all the AF configurations (compare Tables M and |xil). 
Considering additional layers leads to further symmetry 
breaking and renders the distinction between the config- 
urations impossible. The distinction between the cases A 



VI 



and C is possible in confs. a) and k) (compare Tables 
and XII] ) . Besides, in most configurations it is possible 
to decide if these additional layers play any role (compare 
Tables |XI| and |XIIlD. 



D. Distorted oxygen sublattice 

Due to the strong charge-transfer between nickel and oxy- 
gen in NiO the sublattices may be distorted. This effect 
can lower the symmetry of the surface. A point-charge 
model calculation by Iguchi and Nakatsugawa |^l| pre- 
sented a shift of the oxygen sublattice ( "rumpling" ) in the 
direction perpendicular to the surface. Their method did 
not show any in-plane displacement and thus no change 
of the surface symmetry. However, if the "rumpling" 
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also has an in-plane component, i.e. if the oxygen atoms 
are displaced also in the x and y directions, it will also 
have a considerable effect on the symmetry of the crystal 
surface. For this paper, we have chosen a distortion that 
can lower the symmetry of the surface and besides can be 
represented within one conventional unit cell. The config- 
urations considered here are "ferrol" , "ferro2" , "ferro4" , 
a), b), c), e), f), i), k), m), o), p), and r) for the (001) 
surface (see Fig. |l|), "ferrol", ferro3", "ferro5", a), c), f), 
i), and k) for the (111) surface (see Fig. ||), and all con- 
figurations depicted in Fig. ^ for the (110) surface. Other 
depicted spin structures form domains of these configu- 
rations and are not referred to in this subsection nor in 
the tables concerning the current subsection. 

As will be shown later, the best conditions for the de- 
tection of this kind of distortion are presented by the 
(110) surface. The (111) surface could show equally good 
possibilities if only a monolayer of magnetic atoms is 
present. 

In the presence of an oxygen sublattice distortion, the 
chemical unit cell is also doubled. This effectively means 
that magnetic unit-cell-doubling (describing the fact that 
the magnetic unit cell is twice as big as the chemical one) 
is lifted. In general, taking into account distorted oxygen 
atoms in the paramagnetic phase does not lower the sym- 
metry of the problem. The exception is the (111) surface, 
where the six-fold axis is replaced by the three-fold one. 

In the case of the distorted oxygen sublattice, the sym- 
metry group for each configuration is a subgroup of the 
corresponding "non-distorted" configuration, i.e. of the 
corresponding spin configuration in the case where the 
oxygen atoms are not considered. As in case C we dis- 
play only the paramagnetic phase in Fig. || to depict the 
atom positions. All the spin configurations are the same 
as for the corresponding surfaces in case A, and the spins 
are assume d to b e equivalent. 

As Table XIV shows, six different responses can be ex- 
pected from the (001) surface. The paramagnetic surface 
will give a characteristic response. The second group is 
formed by the confs.: a), b), e), o), and "ferrol". Al- 
though confs. a) and o) have elements different from the 
remaining configurations in this group, this fact corre- 
sponds simply to rotating the sample by 90° with respect 
to the z axis. Confs. c) and f) reveal all tensor elements 
without relations between them. Confs. "ferro2", i), k), 
and m) reveal all tensor elements with some relations. 
The only difference between conf. m) and others from 
this group is like for the previous group a 90° rotation 
with respect to the z axis. Another group consists of 
conf. p) alone. It reveals the same tensor elements as 
the paramagnetic phase, but certain relations between 
tensor elements are broken due to a lower symmetry of 
the conf. p). The confs. r) and "ferroS" form the last 
group. All the configurations but k) and "ferroS" can 
be distinguished from those of case A (compare Tabs. || 
and XIV). However only confs. c) and g) can be distin- 



guished from case C (compare Tabs. pC|and pav|). Thus, 
only in these configurations it will be possible to detect 
oxygen sublattice distortions by SHG. 

The SHG response types for the (110) surface are pre- 
sented in Table XV. One can observe that only configu- 



rations c), f) and i) give rise to new (compared to case A, 
Table ^V| ) tensor elements. Compared to case C (Table 
XI), confs. c), f), and i) bring new tensor elements, and, 
surprisingly, confs. a) and g) have less tensor elements, 
due to higher symmetries in the case D. Consequently, 
the confs. a), c), f), g), and i) allow for an unambiguous 
determination of the oxygen sublattice distortion from 
the (110) surface. The possibility to distinguish different 
configurations is rather limited. 

Oxygen sublattice distortion similar to the one pre- 
sented in Fig. H for a (111) surface was found by Renaud 
et al. and calculated by Gillan in M2O3 materials 
(M — Al, Fe). Since the nonmagnetic sublattice symme- 
try group has an infiuence on SHG this distortion can 
be detected also on surfaces of fee crystals. In the previ- 
ous cases A and C we divided the study of (111) surfaces 
in two subcases, considering either one or more atomic 
layers. Taking into account a distorted oxygen sublat- 
tice leads us immediately to the subcase of "more atomic 
layers" . It is caused by the fact that the oxygen and mag- 
netic atoms belong to mutually exclusive planes. The re- 
sulting SHG response types are listed in Table XVl. For 



the AF and ferromagnetic phases, all tensor elements are 
allowed for every configuration. Thus SHG cannot detect 
the magnetic phase of the surface nor distinguish differ- 
ent configurations. Only confs. paramagnetic, "ferroS", 
"ferroS" , and d) allow to decide unambiguously whether 
the oxygen sublattice is distorted or not (compare Tabs. 



VI, XIH, and XVl) 



For both the (001) and (111) surfaces, the symmetry 
groups of case D appear to be the subgroups of the cor- 
responding configurations of case C. This means that the 
oxygen sublattice distortion makes some (one half of all) 
magnetic atoms distinguished as in case C, even though 
we did not apply this distinction explicitly in case D. On 
the other hand, the symmetry groups of the case D dif- 
fer essentially from those of case B. This is caused by 
the difference in distortions assumed in these cases: the 
rhombohedral one in case B and rotation-like in case D. 



E. Domain imaging 

For simplicity, we will consider here only surfaces de- 
scribed hitherto by the case A of our analysis. In this 
case, for AF surfaces, no 180° domains can be expected 
due to the presence of magnetic unit-cell doubling. The 
allowed domains can be detected by surface-sensitive 
SHG under the following two conditions. 

First, domains can be imaged by our method only if 
they manifest themselves at the surface, i.e. if the surface 
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spin ordering changes while passing from one domain to 
another |Q. It is necessary to note, however, that the 
spin orderings for different domains must befong to the 
same configuration in the sense of our classification. We 
do not consider it as a domain structure if one portion of 
the surface is in one configuration and another portion 
is in a different configuration. Under such conditions, 
we can encounter two different types of domains: 90° 
domains (for the (111) surface they are rather 60° do- 
mains), resulting from the rotations around z axis, and 
the second type (called by us mirror-domains, character- 
istic for antiferromagnets), where spins point along the 
same axis in all domains, but the ordering is still different 
(they are no 180° domains!). The tables contain the com- 
plete information about the parity of tensor elements in 
mirror-domain operations, and also for 90° type domains, 
but not for 60° domains. The 90° type domains will be 
addressed later on. In the mirror-domain structure, the 
magnetic point group describing the configuration must 
lack an operation that, while belonging to the (nonmag- 
netic) point group of the system and leaving the spin axes 
invariant, only flips some of the spins. Note, the flipped 
subset of the spins must be antiferromagnetically ordered 
in itself. Configurations, the symmetry groups of which 
lack one of these operations can reveal surface domains, 
related to each other by this operation. 

For an illustration we choose the configuration c) of 
the (001) surface (see Fig. p. The spins point along the 
X axis. Thus operations leaving the axis invariant are 2^, 
2y and 2^ . Of them, 2^ and 2y are absent in the magnetic 
point group of the considered configuration (see Tab. ||, 
conf. c), and Tab. ||). The flipped subset of spins consists 
of the four outer spins for the 2^ operation, and of the 
central spin for 2y (see Fig. ^ b) and c), respectively). 
In fact, there are two domains possible in this configura- 
tion: one with the spins kept invariant under translations 
by the vector (— f , f ,0) (this domain is shown) and the 
other with the spins kept invariant under translations by 
the vector (|, |, 0). Here, a denotes the lattice constant. 
These domains are depicted in Fig. ^. 

The second condition for domain imaging is an inter- 
ference. It can be created internally by different elements 
of the tensor x^'^'^^ ot by external reference [^J46| . The 
interfering elements should be of a similar magnitude for 
the largest possible image contrast. Group theory, how- 
ever cannot account for the amplitudes. With external as 
well as internal reference, a tensor element that changes 
its sign under the reversal of the antiferromagnetic order 
parameter L is necessary. Actually, every L dependence 
of x^^"-* can be represented by splitting the tensor ele- 
ments into odd and even ones in L; even if a tensor ele- 
ment is not purely odd or even we can always decompose 
it according to: 

(2i^) _ (2u),odd (2uj),even , r\ 

i.e. a tensor element consists of parts which are odd and 



even in L, respectively. In a system with many terms 
of that kind the possibility of detecting domains may be 
limited, since they can influence the signal with opposite 
sign, thus diminishing the interference. In highly sym- 
metric structures, such as an fee crystal, the situation is 
more comfortable: every tensor element is either odd or 
even in L (see Appendix). By the appropriate set of ex- 
periments an element can be singled out and give a clear 
image of AF domains. 

As an example we consider tensor elements that are 
present in all the phases, e.g. xizz^: they are even in 
the magnetic order parameters L and M, for the AF and 
ferromagnetic phases, respectively. The tensor element 
XzaJy', present for example in the previously discussed 
conf c) of the (001) surface (see Fig. is odd, since 
it changes its sign under the operation 2^ transforming 
one domain into another. For other configurations other 
tensor elements and operations can be found. In the 
discussed configuration both these elements are present, 
we have intensity contributions proportional to (xizz')^, 

(xi^y^)^ and xizz^ ■ xiiw\ due to the square in Eq. (2). 
As a result, one obtains an interference: 

Ip - ... + (X?zz')^ + ± 2xg%^ ■ XS^ + ... 

(6) 

where stands for one domain, "-" for a different one. 

Now, we turn to the 90° domain structure. Again, 
we take the conf. c) of the (001) surface as an exam- 
ple. The operation connecting the domains is 4z. Under 
this operation, the tensor element x'^^xy changes its sign, 
thus again we have an interference which renders the do- 
main imaging possible. This tensor element is even in 
the domain operation 2xy (which is equivalent to the su- 
perposition of 2x and 4^), which means that domains 
related to each other by this operation cannot be imaged 
using this particular tensor element. Similarly, if a ten- 
sor element is odd in one domain operation and even in 
another, it must be odd in their superposition. Concern- 
ing the 60° domains for (111) surfaces, the parity of the 
tensor elements must be treated more carefully, as indi- 
cated already in . We can still define three "twofold" 
operations, and each of them has its own set of odd and 
even tensor elements. The sets corresponding to differ- 
ent of those operations are not mutually exclusive, i.e. a 
tensor element is usually shared among different parities. 
In this way, this tensor can be positive in one domain, 
negative in the second, and zero in the third one. Thus, 
the existence of a well defined parity of tensor elements is 
necessary for domain imaging, but not sufficient for 60° 
and 120° domain structure. 

This unleashes an interesting question of the antifer- 
romagnetic order parameter. There are as many order 
parameters as different domain structures for a given con- 
figuration. For 60° and 120° domain structures, the AF 
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order parameter must be a vector, while for mirror do- 
mains it is a number. For 90° domains it is can be also a 
number, since there are only two 90° domains. The vec- 
torial order parameter transforms itself under the domain 
operation like a usual vector. 

It is necessary to mention at this point that taking 
into account the spin structure in the second layer would 
not change the validity of the analysis presented in this 
subsection. The only modifications would result from 
addressing bulk domains rather then surface domains, 
and the symmetry of the AF configurations would be 
changed. Yet it would still be possible to find domain 
operations as well as odd and even tensor elements lead- 
ing to interference and AF domain contrast. However the 
possibility to identify each of the domains may be limited 
in some cases due to the increased number of domains. 



IV. POSSIBLE EXPERIMENTAL SETUPS 

In this section, we propose and discuss possible experi- 
mental setups for the detection of AF configuration and 
the imaging of AF domains from low-index surfaces of 
NiO that exhibit magnetic unit-cell doubling in contrast 
to bulk CraOa |,|2^. We propose an experimental setup 
for the detection of antiferromagnetism in the following 
way: both the incident and reflected beams may lie in the 
xz plane (optical plane), and form the angle with the 
z-axis (normal to the sample surface). In the plane per- 
pendicular to the outgoing beam axis, the electric field 
of the second-harmonic generated light has two compo- 
nents, E^^^ and Ef'^'^ , given by the formulae 

|£;(2'^)| = I cos'dEf'^^ - sini?M2'^)| 

p(2t^)| ^ |p(2c^)| ^ ' 

\-^s I \-^y I 

E^^\ and E^z^^ are the components of the elec- 
tric field resulting from SHG in the coordinate system of 
the sample. The dependence of these components on the 
input electric field is indicated by the tensor ^"^^^ . The 
aim of the experiment is the determination of vanishing 
and nonvanishing tensor elements. The easiest way to do 
this is to analyze the output signal intensity as a function 
of the input polarization in both output polarizations s 
and p, for a fixed angle of incidence and reflection. The 
dependence of the output second-harmonic electric field 
on the input polarization is schematically displayed in 
Fig. [To|a)-c) for all tensor elements. The intensity of SHG 
light is the square of the linear combination of these par- 
tial responses. An example of the intensity dependence 
on the input polarization is presented in Fig. [To|d). The 
intensity need not be symmetric with respect Xo ^ — 90° , 
this results from the influence of the electric field depicted 
in Fig. pO|c) . The coefficients of the mentioned combina- 
tion are the products of the x'^^"-' tensor elements and the 



corresponding Fresnel coefficients, according to Eq. (2). 
Thus performing a best fit of these coefficients to the ex- 
perimental results will give (after taking into account the 
Fresnel and geometrical coefficients, known for the given 
experimental geometry and material [^) a set of non- 
vanishing elements of the x*-^"-* tensor. Thus for instance, 
the magnetic phase can be determined. 

Concerning another experimental geometry, with input 
polarization fixed and intensity measured as a function 
of the output polarization, it is possible to determine 
whether the nonlinear Kerr effect takes place. For in- 
stance, with the input polarization ip = 90°, the output 
electric field is given as follows [^ : 

cos<fiA^ie)xgl}B2W + A,ie)xg%^B2m (8) 

As the result, maximum of the intensity is for $ ^ 90°, 
if at least one of the tensor elements x'xyy or x^^m does 
not vanish. Actually, tensor element x^^yy^ is even in all 
the investigated order parameters, but the tensor element 
X^xyy can be odd. For such configurations the Kerr effect 
(change of polarization caused by inversion of the mag- 
netic order parameter) takes place. Thus, it is possible 
to determine which tensor elements are associated with 
the spin-orbit coupling. 

The geometry with p polarization of the reflected SHG 
light seems to be less useful, since there the tensor ele- 
ment xizz"* is always present, regardless of the configura- 
tion. Besides, this polarization mixes the x^x.^"^ and xi^J^"* 
tensor elements. This mixing, however, can be tuned by 
varying the angle of incidence ■& and taking into account 
the influence of the Fresnel coefficients. For smaller {} 
only the Xa^"^ elements are important, while for larger 
the xi^J^'' dominate. If the experiment does not show any 
difference for these two situations, the tensor elements 
must be related. This is the possibility to distinguish 
the configurations with some relations between the ten- 
sor elements from those without such relations. On the 
other hand, the p polarization is useful for AF domain 
imaging. Thus one of the experimental possibilities is to 
carry out the measurements first in s polarized outgoing 
SHG light to make sure that the material is in the AF 
phase and determine its spin configuration. Then a sec- 
ond measurement in p polarization can be performed for 
the domain imaging. 

V. CONCLUSIONS 

Already a short look at the presented tables shows that 
our method works best if the paramagnetic phase is of 
high symmetry, since then a wide variety of different sym- 
metries exists which may be broken by different spin con- 
figurations. In other words, there is enough room for dif- 
ferent new tensor elements to appear along with different 
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spin ordering under these circumstances. In general, this 
is the main reason why only nonlinear optics is suited 
for the detection of antiferromagnetism and the imaging 
of AF domains. The linear susceptibility tensor has too 
low a number of elements for these purposes in order to 
produce unambiguous results. Similarly, among the con- 
sidered surfaces, the (110) surface is the least useful for 
the analysis as it yields ambiguous signal interpretations 
due to its low symmetry in the paramagnetic phase, and, 
on the other hand, very similar symmetries in all the AF 
configurations. 

The (001) and (111) surfaces present alike possibili- 
ties of distinction between the cases. If more than one 
monolayer is involved, however, the (111) surface will 
give the same response in the cases A (all atoms equiv- 
alent) and C (two kinds of magnetic atoms). Both the 
(001) and (111) surfaces also allow for the determination 
of the spin structure, provided the case is known. The 
(111) surface in the case D (oxygen sublattice distortion) 
is an exception - all the AF configurations produce the 
same response. It is possible, however, to determine the 
phase of the material. 

The case D appears to be a subgroup of the case C, 
i.e. all the magnetic point groups describing the configu- 
rations of the case D are subgroups of the corresponding 
ones in the case C. The only exception is the (110) sur- 
face. This inclusion means that the oxygen sublattice 
distortion makes some (one half of all) magnetic atoms 
distinguished as in case C, even though we did not apply 
this distinction explicitly in case D. 

From the fact, that the influence of oxygen sublattice 
distortion (case D) is not detectable in the paramagnetic 
and ferromagnetic phases it follows that only antifer- 
romagnetic ordering can give an extensive information 
about the structure of the surface. It is the magnetic 
atoms and their magnetism which reveal the presence 
and position of oxygen. 

Our short analysis of an AF bilayer structure (surface 
(001)) indicates very similar features to the (GDI) mono- 
layer. There exists a possibility to distinguish AF con- 
figurations from each other, and a certain possibility to 
detect the magnetic phases. Furthermore, introducing 
the second atomic layer does not affect the possibility to 
image AF domains. 

Concerning the magnetic phases, configurations and 
cases considered in this paper, some a priori information 
about the structure is needed in order to draw unique con- 
clusions from the experimental results. For the detection 
of the phase and the spin configuration this additional 
information is the case (A, B, C, D). Vice versa, the case 
(for instance a possible distortion of the oxygen sublat- 
tice) can be determined if one knows the configuration 
(and if it had been previously deduced that the investi- 
gated material is antiferromagnetic). Actually, in most 
measurements of AF spin structures some a priori knowl- 
edge is required. For example, in experiments by Fiebig 



et. al. |2l| such a prerequisite is the assumption of the 
AF spin-flop phase of the material. In both experimen- 
tal approaches mentioned here the (001) surface seems to 
provide the best possibilities of drawing valuable conclu- 
sions, while the (110) surface is the least suitable in that 
respect. 

Finally, our paper demonstrates that the AF domain 
imaging is possible even in the presence of magnetic unit- 
cell doubhng. Thus optical SHG, unlike linear optics, is 
able to image AF surface domains. For most AF conflg- 
urations, there are more than one surface domain struc- 
tures. The rule stating that the number of domains is 
equal to the number of symmetry operations in the para- 
magnetic phase divided by the number of symmetry op- 
erations in the magnetic phase is applicable also for an- 
tiferromagnets (thus, with unit-cell doubling the number 
of domains is reduced by a half). However, not all the 
domains can be imaged at the same time. 



APPENDIX 

ON THE GROUP-THEORETICAL ANALYSIS OF 
MAGNETIC SYSTEMS 

In this Appendix, we would like to address some par- 
ticularities of our group-theoretical analysis. The first 
general remark is that although symmetry analysis can 
provide us with a set of nonvanishing tensor elements for 
a given configuration, but cannot give any information 
about their magnitude. This equally applies to the dis- 
tortion effects, as treated e.g. in Ref. p9| . 

Another interesting issue is the behavior of the ten- 
sor elements with respect to the AF order parameter L 
(for ferromagnetic phases L should be replaced by the 
magnetization M), i.e. the parity of tensor elements. In 
general, a tensor element consists of even and odd parts 
with respect to L, as shown in Eq. (^). In systems with 
high symmetry, it is possible to describe an operation 
which reverses L (or M) by a spatial operation I. The 
operation I belongs to the point group of the system, but 
not to its magnetic point group. The application of this 
operation to a tensor element will change its sign (keep 
it invariant) if this element is odd (even) in L. Conse- 
quently, each tensor element can be either odd or even 
in L, a mixed behavior is forbidden. Actually, the parity 
of a given tensor element is a function of the chosen op- 
eration I. In most antiferromagnetic configurations more 
than one operation leading to different domain structures 
are possible (this means, more than one order parame- 
ter can be defined). For example, for (001) surface one 
has 4^ rotations leading to different domains in addition 
to the eventual mirror-structure. For the (111) surface, 
there are three domains resulting from the rotations with 
respect to the z axis alone. For some configurations, they 
exist in addition to the mirror-domains. 
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This whole analysis of the parity of the tensor elements 
cannot be performed for the systems with a lower sym- 
metry, where it is impossible to find an operation I de- 
scribing the inversion of L or M, a mixed behavior is then 
allowed. Note, the presence of dissipation (redistribution 
of response frequencies) does not influence the above con- 
sideration. In general, dissipation in frequency space is 
responsible for the mixing of the real and imaginary parts 
in the tensor elements, while point-group symmetry gov- 
erns the (non)existence of tensor elements purely odd or 
even in the magnetic order parameters L or M. 
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FIG. 1. Spin configurations of an fee (001) surfaee. Except for eonfs. "ferro4" and o) - r), the arrows always indicate in-plane 
directions of the spins. In confs. "ferrol" and o) - r) ((g)) denote spins pointing along the positive (negative) z-direction, 
respectively. 




FIG. 2. Top view of a spin structure on a (001) surface. The dashed line depicts a conventional unit cell, while the solid 
one outlines the primitive unit cell. 
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FIG. 3. Spin structure of an antiferromagnetic (001) bilayer constructed from a shift of the monolayer along the positive x 
(y) axis. Filled (empty) circles represent the topmost (second) layer. On the right hand side the conventional unit cells for the 
resulting bilayer structure are presented. Here, conf. a) of the (001) monolayer serves as an example. 
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FIG. 4. Spin configurations of an fee (110) surface. Except for confs. "ferro3", g), h), and i), the arrows always indicate 
in-plane directions of the spins. In confs. "ferroS", g), h), and i) ((g)) denote spins pointing along the positive (negative) 
z-direction, respectively. 
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FIG. 6. Structure of the (001) and (111) surfaces of a fee crystal with a rhombohedral distortion in the paramagnetic phase. 
Note the changed orientation of the coordinate system for the (001) surface. 
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FIG. 7. Surface structure of the non-equivalent magnetic atoms case in the paramagnetic phase. Pictures present the (001), 
(110), and (111) surfaces, respectively. Filled and empty circles represent the two kinds of magnetic atoms. Note, the fragment 
representing the (111) surface does not show the conventional unit cell but a bigger set of atoms in order to give a clear idea 
about the surface structure. 
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FIG. 8. Surface structures of the case with a distorted oxygen sublattice (white circles). Pictures present the paramagnetic 
phase of (001), (110), and (111) surfaces, respectively. Note, the fragment representing the (111) surface does not show the 
conventional unit cell but a bigger set of atoms in order to give a clear idea about the surface structure. 




FIG. 9. Two surface mirror domains for an AF configuration - panels b) and c) depict the same AF domain, related to the 
panel a) by different mirror operations. 
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TABLE I. Details of SHG response types. We denote Xijk^ by ijk. Odd elements are in bold if a domain operation exists. 



key point group symmetry operations 



domain operation non-vanishing tensor elements 



a 
b 



4mm 
m 



mm2 
2 



mm2 



mm2 
m 



1 O A.A o o o o 
^2: -1-^2? ^yj ^xyi ^—xy 



1,2^ 



1,2,, ±4, 



1,2, 
1,2„ 



1,2- 



1, 22 , 2a; , 2y 

1 9 



1,2, 
1,2, 



2z) 2y 

4 9 

22 , 2_a;y 



4z , 2^, 



9 9 



±4, 



9 9 



2x; 2y 
i42: , 2xy ; 2_; 
i4z , 2a; , 2y 
2z 5 2a; 

42 , ^xy 

2, 



2x 



±4, 



4z , 2y 



2z , 2ii 



23; , 2i( 



22 , 



xxz = xzx = yyz = yzy, zxx = zyy, zzz 

xzx = xxz, xxy = xyx, yxx, yyy, yzz, 

yyz — yzy, zxx, zyy, zzz, zyz — zzy 

no information about the parity 

XXX = -yyy, xyy = -yxx, xzz = -yzz, 

xyz = yxz = xzy = yzx, xxz = xzx = yyz = yzy, 

xxy = -yyx = xyx = -yxy, zxx = zyy, zzz, 

zxz = zzx = -zyz = -zzy, zxy = zyx 

XXX = -yyy, xyy — yxx, xzz — -yzz, 

xyz = xzy = yxz = yzx, xxz = xzx = yyz = yzy, 

xxy = -yyx = xyx = -yxy, zxx = zyy, zzz, 

zxz — zzx = zyz — zzy, zxy = zyx 

xyz = xzy = -yxz = -yzx, 

xzx = xxz = yzy = yyz, zxx = zyy, zzz 

xxz = xzx, yyz = yzy, zxx, zyy, zzz 

xyz = xzy, xxz = xzx, yyz = yzy, yzx = yxz, 

zxx, zyy, zzz, zxy — zyx 

xyz = xzy, xxz = xzx, yyz = yzy, yzx = yxz, 
zxx, zyy, zzz, zxy = zyx 

xxz = xzx = yyz = yz^j, xzy = xyz = yzx = yxz, 

zxx = zyy, zzz, zxy = zyx 

XXX, xyy, xzz, xxz = xzx, yyz = yzy, 

yyx = yxy, zxx, zzz, zzx = zxz 

XXX, xyy, xzz, xxz = xzx, yyz = yzy, 

yyx = yxy, zxx, zzz, zzx = zxz 

All the elements are allowed: 

XXX, xyy, xzz, xyz — xzy, xzx = xxz, 

xxy = xyx, yxx, yyy, yzz, yyz = yzy, 

yzx = yxz, yxy = yyx, zxx, zyy, zzz, 

zyz = zzy, zzx = zxz, zxy = zyx 

XXX, xyy, xzz, xyz = xzy, xzx = xxz, 

xxy = xyx, yxx, yyy, yzz, yyz = yzy, 

yzx = yxz, yxy = yyx, zxx, zyy, zzz, 

zyz = zzy, zzx = zxz, zxy = zyx 

no information about the parity 

XXX = yyy, xyy = yxx, xzz = yzz, 

xyz = yxz = xzy = yzx, xxz = xzx = yyz = yzy, 

xxy = yyx = xyx = yxy, zxx = zyy, zzz, 

zxz = zzx = zyz = zzy, zxy = zyx 

XXX =yyy, xyy =yxx, xzz — yzz, 

xyz = yxz = xzy = yzx, xxz = xzx = yyz = yzy, 
xocy = xyx = yyx = yxy, zxx = zyy, zzz, 
zxz = zzx — zyz = zzy, zxy — zyx 

xxz = xzx, yyz = yzy, zxx, zyy, zzz 
xzx = xxz, xxy = xyx, yxx, yyy, yzz, 
yyz = yzy, zxx, zyy, zzz, zyz = zzy 
All the elements are allowed: 

XXX, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
xxx, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy — zyx 
xyz = xzy, xxz = xzx, yyz = yzy, yzx = yxz, 
zxx, zyy, zzz, zxy = zyx 
xxx, xyy, xzz, xxz = xzx, yyz = yzy, 
yyx = yxy, zxx, zyy, zzz, zzx = zxz 
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p 6mm l,2^,±3^,±6^,6(2x) 

q 6 1,2,, ±3,, ±6, 2,, 2, 

r 3m l,±3z,2y,2s(xy),'2s(-xy) 

S 1 1 2y 

t m l,2j; 

u 3 1,±32 2y 

w 1 1 



TABLE II. SHG response for all spin configurations of the 
(001) surface of a fee lattice For the detailed descrip- 

tion of the response types see Tab. ^ The configurations are 
depicted in Fig. |l|. 

configuration key (response type) 



para a 

ferrol b 

ferro2 c 

ferrol d 

AF: 

a), b), e), o) e 

c), f) f 

i), k), m), p) g 

r) a 



TABLE III. SHG response for all spin configurations of the 
(001) surface of a fee lattice, with the spin structure of the 
second layer taken into account. For the detailed description 
of the response types see Tab. ^ For the confs. see Fig. ^ 

configuration key (response type) 



para a 

ferrol b 

ferro2 c 

ferrol d 

AF: 

ax), ox) h 

ay), oy), r) e 

bx), by), ex), ey) b 

c), fx), fy) i 

i) j 

k) f 

m), p) c 



xxz = xzx — yyz — yzy, zxx — zyy, zzz 

xyz = xzy — -yxz — -yzx, xxz — xzx — yyz — yzy, 

zxx — zyy, zzz 

zxx — zyy, xxz — xzx — yyz — yzy, zzz, 
XXX = -xyy = -yxy = -yyx 
All the elements are allowed: 

XXX, xyy, xzz, xyz = xzy, xzx = xxz, xxy = xyx, 
yxx, yyy, yzz, yyz = yzy, yzx = yxz, yxy = yyx, 
zxx, zyy, zzz, zyz = zzy, zzx = zxz, zxy = zyx 
XXX, xyy, xzz, xxz = xzx, yyz = yzy, 
yyx = yxy, zxx, zyy, zzz, zzx = zxz 

XXX = -xyy = -yxy = -yyx, xyz = xzy = -yxz = -yzx, 
xzx = xxz = yyz = yzy, xxy = xyx = yxx = -yyy, 
zxx = zyy, zzz 
All the elements are allowed 



TABLE IV. SHG response for all spin configurations of the 


(110) surface of a fee lattice llq 


. For the detailed description 


of the response types see Tab. 


n. The configurations are de- 


picted in Fig. ^. 




configuration 


key (response type) 


para 


k 


ferrol 


I 


ferro2 


m 


ferro3 


n 


ferrol 


o 


AF: 




a), b), c), g) - 1) 


k 


d), e), f) 


n 


TABLE V. SHG response for all spin configurations of the 


(111) surface of a fee lattice 


Iq]. Only one monolayer is 


taken into account. For the detailed description of the re- 


sponse types see Tab. |. The configurations are depicted in 


Fig. 1 




configuration 


key (response type) 


para 


P 


ferrol 


1 


ferro3 


o 


ferroS 


q 


AF: 




a), i), k) 


k 


c),f) 


n 
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TABLE VI. SHG response for all spin configurations of the 
(111) surface of a fee lattice [^. More monolayers are taken 
into account. For the detailed description of the response 
types see Tab. |. The configurations are depicted in Fig. 

configuration key (response type) 



para r 

ferrol s 

ferro3 t 

ferroS u 

AF: 

a), i), k) t 

c), f) u 



TABLE IX. SHG response for all spin configurations of 
the (111) surface of a fee lattice, distorted to a rhombohedral 
structure. More monolayers are taken into account. For the 
detailed description of of the response types see Tab. |. For 
the surface structure see Fig. ^ for the spin configurations 
see Fig. |. 

configuration key (response type) 



para t 

ferrol, ferro2, ferro4, ferroS s 

ferroS t 
AF: 

a), i), k) s 

b)-h),j), l),m) t 



TABLE VII. SHG response for all spin configurations of 
the (001) surface of a fee lattice, distorted to a rhombohedral 
structure. For the detailed description of the response types 
see Tab. 0. For the surface structure see Fig. ^, for the spin 
configurations see Fig. [l| 



configuration 



key (response type) 



para 
ferrol 
ferro2 
ferro3 
ferro4 
AF: 
b) - h), o) 
■ n), p) - r) 



k 
m 
o 
1 
n 

n 
k 



TABLE VIII. SHG response for all spin configurations of 
the (111) surface of a fee lattice, distorted to a rhombohedral 
structure. Only one monolayer is taken into account. For the 
detailed description of the response types see Tab. ^. For the 
surface structure see Fig. ^ for the spin configurations see 
Fig. i. 



configuration 

para 
ferrol, ferro4 
ferro2 
ferroS 
ferroS 
AF: 
a),k) 
b)-,i), 1), m) 



key (response type) 



k 
1 

m 
o 
n 

k 
n 



TABLE X. SHG response for all spin configurations of the 
(001) surface of a fee lattice, with one atom distinguished. 
For the detailed description of the response types see Tab. ^ 
For the surface arrangement see Fig. |^ For the confs. see 
Fig. g 



configuration 


key (response type) 


para 


a 


ferrol 


b 


ferro2 


c 


ferro4 


d 


AF: 




a), o) 


h 


b),e) 


b 


c) 


f 


f) 


i 


i), m), p) 


e 


k) 


j 


r) 


d 


TABLE XL SHG response for all spin configurations of the 


(110) surface of a fee lattice, with one 


atom distinguished. For 


detailed description of response types 


1 see Tab. ^ For the sur- 


face arrangement see Fig. Q For the 


confs. see Fig. ^. 


configuration 


key (response type) 


para 


k 


ferrol 


1 


ferro2 


m 


ferroS 


n 


ferro4 


o 


AF: 




a) 


1 


b),c), h),i),k),l) 


k 


d) 


Im 


e), f), g) 


n 


■i) 


o 
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TABLE XII. SHG response for all spin configurations of 
the (111) surface of a fee lattice, with one atom distinguished. 
Only one monolayer taken into account. For the detailed de- 
scription of the response types see Tab. |. For the surface 
arrangement see Fig. |^ For the confs. see Fig. 

configuration key (response type) 



TABLE XIV. SHG response for all spin configurations of 
the (001) surface of a fee lattice, with a distortion of oxygen 
sublattice. For the detailed description of the response types 
see Tab. ^ For the surface arrangement see Fig. ^ For the 
confs. see Fig. 0. 

configuration key (response type) 



para 


P 


ferrol 


1 


ferroS 


o 


ferro5 


q 


AF: 




a) 


1 


c),f) 


m 


i) 


o 


k) 


n 



TABLE XIII. SHG response for all spin configurations of 
the (111) surface of a fee lattice, with one atom distinguished. 
More monolayers are taken into account. For the detailed de- 
scription of the response types see Tab. |. For the surface 
arrangement see Fig. |^ For the confs. see Fig. 



configuration 



key (response type) 



para 
ferrol 
ferroS 
ferro5 
AF: 
, c), f), k) 
i) 



3r 
s 
t 
u 

s 

t 



para 




a 


ferrol 




b 


ferro2 




c 


ferro4 




u 


XV- 

l\r . 






a), o) 




h 


b),e) 




b 


c),f) 




i 


i),k) 




c 


m) 




j 


P) 




e 


r) 




d 


TABLE XV. SHG 


response for all spin 


configurations of 


the (110) surface of a 


fee lattice, with oxyg 


;en sublattice dis- 


torted. For the detailed description of the response types see 


Tab. |. For the surface arrangement see Fig 


. |8[ For the confs. 


see Fig. y. 






configuration 


key (response type) 


para 




k 


ferrol 




1 


ferro2 




m 


ferroS 




n 


ferro4 




o 


AF: 






a), b), g), h), k), 1) 




k 


c) 




o 


d), e), i), j) 




n 


f) 




m 



TABLE XVI. SHG response for all spin configurations of 
the (111) surface of a fee lattice, with oxygen sublattice dis- 
torted. For the detailed description of the response types see 
Tab. Q. For the surface arrangement see Fig. ^ For the confs. 
see Fig. |. 

configuration key (response type) 



para u 

ferrol, ferroS w 

ferroS u 
AF: 

All confs. w 
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